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Abstract 

Recent experiments on the quasi-lD antiferromagnets R2BaNiO^ 
( R = rare earth ) have shown the existence of purely ID Haldane gap 
excitations propagating on the Ni chains. Below an ordering temper- 
ature, the gap excitations survive and coexist with the conventional 
spin waves in the ordered phase. We construct a model mixed-spin 
system in 2D for which the ground state can be exactly specified. Us- 
ing the Matrix Product Method, we show the existence of Haldane gap 
excitations in the ordered phase. We consider different cases of order- 
ing to study the effect of ordering on the degeneracy of the Haldane 
gap excitations. 
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I. Introduction 

Haldane^ in 1983 predicted that a Heisenberg antiferromagnetic (AFM) chain 
of integer spins has a gap in the excitation spectrum. Later, Haldane's con- 
jecture was verified in several theoretical as well as experimental studies. 
In a quasi- ID magnetic system, the integer spin chains are coupled via weak 
exchange interactions. If the strength of the coupling is very small, there 
is no possibility of ordering and individual spin chains remain in the disor- 
dered Haldane gap phase at all temperatures. In materials like CsNiCl^, 
the inter-chain coupling is sufficiently strong to give rise to an AFM ordered 
phase below the Neel temperature T^-^'^ For temperature T > Tj^, the three 
Haldane gap modes can be observed in experiments but as T approaches Tjv, 
the gap vanishes at the 3D magnetic zone centre. In the ordered phase, two 
of the three modes become conventional gapless spin waves, while the third 
'longitudinal' mode develops a gap. 

For quasi- ID Haldane gap systems, there now exists a third possibility. 
Recent inelastic neutron-scattering experiments^"^ on the mixed-spin anti- 
ferromagnets R2BaNiO^ {R — Pr, Nd or NdxYi^^), provide clear evidence 
that the Haldane gap excitations propagating on the Ni chains survive even 
in the ordered phase. Unlike in systems like CsNiCl^, the Haldane gap does 
not become zero as T approaches T^. In the ordered phase, the gapped Hal- 
dane modes coexist with the conventional spin waves characteristic of the 
ordered phase.The Haldane gap excitations maintain their ID character in 
the entire phase diagram. In this paper, we construct a model mixed spin 
system in 2D for which the ground state can be exactly specified. AFM long- 
range order exists in the ground state with gapless spin waves as excitations. 
The system also consists of integer spin chains, the ground states of which 
are the Valence Bond Sohd (VBS) states. Haldane gap excitations can be 
created in the chains in the ordered phase of the mixed-spin system. We also 
discuss how the gap changes as the nature of the ordering changes. 

II. Ground state and excitations 

Our model mixed spin system has an underlying lattice structure shown in 
Fig. 1. It basically consists of a square lattice at each site of which a spin- 
1/2 (A-spin, sohd circle) is sitting. The centre of each square plaquette is 
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occupied by a S=l spin (B-spin, solid square). The B-spins along a chain 
in the horizontal direction interact via nearest-neighbour (n.n.) interactions 
but the individual B-spin chains do not interact with each other. The B- 
spins interact with the n.n. A-spins along the diagonal lines in Fig. 1. The 
A-spins interact with n.n. A-spins along the horizontal and vertical bonds of 
the square lattice. In R2BaNiO^ also, the S=l Ni chains arc decoupled and 
the sites with half-odd integer spins are positioned between the chains. 
The ^23" spins form a network via Ni-O-R and R-O-R superexchange routes. 

We now demand that the B-spin chains arc in Valence Bond Solid (VBS) 
ground states. Affleck, Kenedy, Lieb and Tasaki (AKLT)^ constructed a 
Hamiltonian for which the VBS state is the exact ground state of an integer 
spin chain. We follow the same principle to construct the Hamiltonian of 
our model. Fig. 2 shows two successive rows of A and B-spins. Consider 
the B-spins to be in a VBS state. Each spin-1 can be considered to be a 
symmetric combination of two spin-l/2's. In the VBS state, each spin-1/2 
forms a singlet with a spin-1/2 at a neighbouring site. This process has to 
be symmetrized so that a spin-1 is restored at each site. In the ground state, 
a valence bond (VB) or singlet covers every link along the integer-spin chain 
so that the total spin of the state is zero. The Hamiltonian for which the 
VBS state is the exact ground state is given by 

H\j, = Y,P^/2{Sa + S, + S,) (1) 

A 

P5/2 is the projection operator onto spin 5/2 for the triangle of spins. The 
sum is over successive triangles of spins with bases along the integer spin 
chain. The spins at sites Bi and B2 have magnitude 1 and the A-spin has 
magnitude 1/2. In the VBS state, the sum of the three spins can never be 
5/2 as a VB exists between the basal spins. Thus, the projection operator 
P5/2 operating on the triangle of spins in the VBS state gives zero. The VBS 
state is the exact ground state of H\q with eigenvalue zero. It can be easily 
verified that the A-spins are free when the Hamiltonian in (1) operates on the 
VBS state of the integer spin chain. The Hamiltonian H\q, on expanding, 
is given by 

H\b = 7^[2 + ^{Sa.S, + + S,.Sa) + 2(51.^5 + S,.S, + S^.Saf] (2) 
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When we consider the full 2D spin system, the Hamiltonian is given by 

Hab = X] ^AB 

i 

The sum is over individual B-spin chains. We now consider the interactions 
between the A-spins. The Hamiltonian Ha is given by the usual Heisenberg 
AFM Hamiltonian defined on the square lattice: 

Ha^JYI Si.Sj (3) 

<ij> 

The spins have magnitude 1/2 and < ij > are n.n.s. The ground state of 
(3), though not exactly known, has long range Neel type of order. The total 
Hamiltonian H is 

H^Ha + Hab (4) 

The ground state of H has the following structure: the B-spin chains are 
in VBS states. The square network of A-spins is in the state which is the 
ground state of the S=l/2, n.n. Heisenberg AFM Hamiltonian on a square 
lattice. The decoupling of the S=l spin and the S=l/2 spin ground states is 
possible because the A-spins are free when Has operates on the VBS states. 

We now consider the excitations in the system. Spin wave excitations 
created in the A network of spins have no cfi'cct on the VBS ground states 
of the chains of B-spins. The spin wave excitations are gapless. Consider 
now the case when the network of A-spins is in the ground state and an 
excitation is created in a B-spin chain. Again, for convenience, we consider 
the two-chain strip of Fig. 2. The A-spins in the top chain are in the Neel 
state Itititi ■ ■ ■) indicative of long range AFM order. We determine the 
excitation spectrum by the method of Matrix Products^"^^ The VBS ground 
state of the chain of B-spins can be written as 

I V'g) = Tr{gi ® 92® • • ■ ® Ql) (5) 

where L is the number of sites (spins) in the chain. The 2x2 matrix gi is 
given by 

''=[V2\l) |0) J 
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The states | 1), | 0), | I) correspond to S"^ = 1, and —1 respectively. One 
can easily verify that Hab \ V'g) = 0. An excited state of the chain of B-spins 
is constructed as [3] 

\ rsik)) ^ ^T.e-''' \ rsj) (7) 

I ip%j) is obtained by replacing gj in the matrix product of the VBS state by 
the matrix Qj, where a = 1, 0, — 1 are the values of the z-component of the 
total spin, 5"^, of the excited state. The matrices have the form 

._/x/2|l) 0\ -|0) V2\l)\ ,_(0 -|0) \ 

(8) 

I il^^Bj) = Trgi ® 52 • • • Qj-i ® gj+i ® ■ ■ ■ ® gL (9) 



^ 1=1 p=l 



We now consider the case a = 1. 



j.i \ j,i \ ^ V ■■■ n^^ ■ ■ ■ 

{na,ma} 

9mxm29m-2,m-i ' ' ' fi'm^mp+i ' ' ' Qmi^mx (H) 

Define three 4x4 matrices G, G and G!' as 

^ ~ ^inimi){ni+-imi+i) — 9nini+i9rnimi+i (12) 

/iiA(2 (npTOp)(np+imp+i) 9npnp+\9mpmp+i 
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The ordering of multi-indices is given by 



/x = 1,2,3,4^ — > 



(11), (12), (21), (22) 



G(nimi){n2m2)Cr(n2m2){n3m3) 

g' c 

(n; m; ) (ni+ 1 mj+ 1 ) (ni+i m(+ 1 ) (n;+2m(+2 ) 

r" r 
TrG^-^G'GP-^+'^G"G^-P 



The matrices G, G' and G" are 
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The eigenvalues and eigenvectors of G are : 

Ai = 3, X2 — A3 

, , ^ 
V2 

V 1 / 

/ \ 

1 



V / 



les) 



|e4) 



A4 = -1 




V 1 / 

/ \ 



1 

V / 



6 



From (14) and in the thermodynamic hmit L— > oo, 



Therefore, 



n=l 
1 



G 



-i+i 



1, 



-3^(--)^'-' 



n 



G 



ei> Ai 



L-p 



I Kik) 



ik(l—p) 



L 

= 3^ 



L L 

/=i p=i 

2 

5 + 3 cos /c 



1 



\p-i\ 



(17) 



The chain of A-spins is in the Neel state | ipNeei)- The total wave function 
I 'i/') of the two- chain system is a product of | V'ij(^)) I '0Areei)-The norm 
{tp I is given by the same expression as in (17). The excitation energy 
measured w.r.t. the ground state energy is 



ii^ \Hab 




(V; 





(18) 



The denominator has aheady been calculated. We now calculate the numer- 
ator. Hab for single chain is given by Eq. (2). The two nearest-neighbour 
sites of the lower chain in Fig. 2 and the intermediate site in the upper chain 
form the vertices of a triangle. If one of the sites of the lower chain has the 
'defect' matrix assosiated with it, then there are four possible states of 
the triangle : 



I "Pij) 
I "P^j) 

One can verify that 



T 

9j ® 9j+i 



i 

9j ® 9j+i 



T 

93 ® 9]+i 
i 

9j ® 9j+i 



(19) 



H 



H. 
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rr I,, X^tll^ -8v^(|10)+|01)) -40 I 11) \ 

^3,j+i\^ij/ ' 20 1^8 I 00) + 4 I 11) + 4 I 11) 8x/2(| 01)+ I 10)) y 

^ 20 V 4^2(1 01)+ I 10)) 8 I 11) J 

1 /-8|00)-4|11)-4|11) -8^2(1 10)+ I 01)) \ 
I 20 1^ 4^2(1 iO)+ I 01)) 8 I 00) + 4 I 11) + 4 I 11) y 

]_( -4v^(|10)+|01)) -8 I 11) A 

^ 20 ^ 8 I 00) +4 I 11) +4 I 11) 4^2(1 01)+ I 10)) y ^ ' 

The up and down spins outside the brackets represent the states of the A- 
spin. 

We represent | ip) by the state 

in which the upper and lower rows contain the states of the upper and lower 
chains respectively. One can verify that 



= 5jyTrG^-^Z{H)G^-^-^ (22) 
where the matrix G is given in (15) and 

(nimi)(n2m2) 

mim2 (j odd) 

{j even) (23) 
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From Eqns. (20) 
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oo , Eq. (22) reduces to 

\ ^Bj( 

= (ei G^-^Z{H) ei 



\{HAB)jj+l\ 



Xr^e,\Zmei) 



= — X 3^-"^ X 60 for j = odd 
20 

= — X 3^~^ X 20 for 7 = even 
20 



Ai and | ei) are given in Eq. (16). From (7), 



{i^\HAB\i^) 
1 



^ \^ \^ „ik(j' -j) c I WNeel \/ tt \ I "^Neel 

1 L , 

= — X -60 + -20 

20 L ^2 2 ' 



3^-2 X 2 



From (18), (17) and (26), one gets 
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uji{k) = -(5 + 3 cos A;) 
9 
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The excitation energy (27) has been determined for a two-chain strip. For 
the 2D spin system, the chain of B-spins is connected to two chains of A- 
spins and the excitation energy is twice of that in (27). The excitation energy 
uj^i{k) for = —1 is the same as Ui{k). The excitation energy a;o(^) for 
5*^ = is 2L0i{k). Thus the gaps in the excitation spectrum for the 2D spin 
system are : 

A.i = ^ , Ao = ^ (28) 

We now consider the case when excitations are created in both the A and 
B-chains. The excitation in the A-chain travels in the network of A-spins as 
a 2D spin wave. The excitation spectrum of the spin wave is gapless. The 
Haldane gap excitation is confined to the B-spin chain. For a ferromagnetic 
ordering of A-spins, one can verify that A+i A_i, i.e., the degeneracy of 
the triplet Haldane modes is fully hfted. 

III. Conclusion 

We have constructed a model mixed spin system for which the ground state 
can be exactly specified. The coexistence of AFM long range order and 
Haldane gap excitations has been explicitly shown. In the ordered phase, the 
degeneracy of the three Haldane gap modes is partially lifted. The transverse 
modes (S^ = ±1) arc still degenerate whereas the longitudinal mode (S^ ~ 0) 
has an excitation gap twice that of the transverse modes. These results are in 
agreement with the experimental observations on the R2BaNiOr, systems. ''"^ 
In these systems, coexistence of gapless spin waves in an ordered phase and 
the Haldane gap excitations confined to Ni chains, has been observed. The 
three Haldane gap modes in an isolated chain of integer spins arc degenerate. 
The lifting of the degeneracy seen in R2BaNiO^ is due to the ordering of the 
rare-earth spins. ^ It has been suggested^^ that in R2BaNiO^ the Ni-spins 
participate in the long range order. This is, however, not so in our model, as 
in the ground state there is a complete decoupling of the A and B-spins. An 
important feature of the R2BaNiO^ systems is the increase in the gap energy 
as the temperature is lowered, i.e., as the system becomes more ordered^^. 
Finite temperature calculations on our model are in progress and the results 
will be reported elsewhere. 



10 



Acknowledgement 

One of the Authors (EC) is supported by the Council of Scientific and In- 
dustrial Research, India under sanction No. 9/15(186)/97-EMR-I. 



11 



References 



ip.D.M. Haldane, Phys. Lett. 93 A, 464 (1983), Phys. Rev. Lett. 50, 1153 
(1983) 

^L. P. Regnault, L Zaliznyak, J. P. Renard and C. Vettier, Phys. Rev. B 50, 

9174 (1994) and references therein 
^G. Path and J. Solyom, J. Phys.: Condens. Matter 5, 8983 (1993) and 

references therein 

^L Affleck, T. Kennedy, E. H. Lieb and H. Tasaki, Phys. Rev. Lett. 59, 799 

(1987); Commun. Math. Phys. 115, 477 (1988) 
^W. J. L. Buyers et ah, Phys. Rev. Lett. 56, 371 (1986); R. M. Morra, 

W. J. L Buyers, R. L. Armstrong and K. Hirakawa, Phys. Rev. B 38, 

543 (1988) 

^K. Kakurai, M. Steiner, R. Pynn and J. K. Kjems, J. Phys.: Condens. 

Matter 3, 715 (1991) 
''A. Zheludev, J. M. Tranquada, T. Vogt and D. J. Buttrey, Phys. Rev. B 

54, 6437 (1996); Phys. Rev. B 54, 7210 (1996) 
^S. Raymond, T. Yokoo, A. Zheludev, S. E. Nagler, A. Wildes and J. 

Akimitsu, |cond- mat /981 1040 
^A. Kliimper, A. Schadschneider and J. Zittartz, J. phys. A 24, L955 (1991); 

Z. Phys.B 87, 281 (1992); Europhys. Lett. 24, 293 (1993) 
^°H. Niggemann and J. Zittartz, Z.Phys. B 101, 289 (1996) 
^^K. Totsuka a nd M. Suzuki, J.Ph ys.: Condens. Matter 7, 1639 (1995) 
^^A. Zheludev, |cond-mat/9707167| and references therein. 



12 



Figure Captions 

Fig. 1. Model mixed spin system consisting of a square lattice of S = 1/2 
spins (A-spins, solid circle). The centre of each square plaquette is 
occupied by a S = 1 spin (B-spin, solid square) . 

Fig. 2. A chain of A-spins (upper row) connected to a chain of B-spins (lower 
row). 
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